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The build-up of exciton populations in resonantly laser excited semiconductors is studied micro- 
scopically. For excitation at the 2s-exciton resonance, it is shown that polarization with a strict 
s-type radial symmetry can be efficiently converted into an incoherent p-type population. As a con- 
sequence, inversion between the 2p and Is exciton states can be obtained leading to the appearance 
of significant terahertz gain. 
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PACS numbers: 71.35.-y, 42.60.Lh, 78.30.-j 



Resonant laser excitation of semiconductors induces 
a coherent interband polarization between conduction- 
band electron and valence-band hole states. Through 
interaction and scattering processes this optical polar- 
ization may be converted into incoherent populations of 
unbound or bound electron-hole pairs (excitons). Even 
though excitonic features of the coherent polarization are 
well understood 0, 0] , the study of its decay into inco- 
herent many-body states is an area of active research. 

Several recent experiments have applied terahertz 
(THz) fields to directly probe the optically generated 
many-body system. So far, this approach has been used 
to detect and monitor conductivity 0, plasmons Q, 
and bound exciton formation At the same time, 

tunable terahertz (THz) sources evolve rapidly ranging 
from free-electron |(| and quantum-cascade lasers |j| to 
sources with difference- frequency generation ||. Since 
the semiconductor quasi-particle excitations strongly in- 
teract with THz radiation, it is an interesting question to 
see if an excitonic system could actually be used to gener- 
ate THz radiation, or even to provide THz amplification, 
i.e. THz gain. 

In order to explore the THz properties of resonantly 
excited semiconductors, we study the build-up of exci- 
ton populations in different quantum states. Besides the 
conversion of excitonic polarizations at the Is or 2s reso- 
nances into incoherent s-type populations, we show that 
Coulomb induced scattering can efficiently convert exci- 
tonic coherences with a strict s-type radial symmetry into 
an incoherent p-type population. We show that this pro- 
cess may even lead to a population inversion between the 
2p and Is states, giving rise to THz probe gain. 

As a model system, we analyze quantum-wire struc- 
tures but we show that the main results are equally valid 
for quantum-well systems. The electronic excitations are 
described by Fermion operators cic^k and oL„) k related 
to an electron with carrier momentum k in the con- 
duction (valence) band. We include the carrier-carrier 
Coulomb interaction as well as the couplings to light 
fields and phonons Q • The dynamic evolution of single- 
particle quantities related to the microscopic polariza- 



tion P k = (<4 ik a c ,k), electron / k = (a^ k a Cjk ), and hole 
/ k = (a„ ik ol k ) densities obey the semiconductor Bloch 
equations [3 
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and similarly for / k . Here, Vic is the Coulomb- 
matrix element whereas e k and f2 k denote the renor- 
malized kinetic energy and Rabi frequency, respec- 
tively. The coupling to the THz field follows from 
the jA term which contains the vector potential Ae a 
with direction e a and the current-matrix element j k = 
e(^- + ^-)7ik • e a with the electron (hole) mass 
m e(h)- The true two-particle correlations stem from 



and T> 



A,// 



q,k',k A / T T \ 

C \,v,l>,\> = A K,k<k''V,k'+q,aA',k-q), — - k>q 

J2 Pz £p*,qA((-Dp*,q + £>jL q ) a A,k a ^k- q ) where D and 
D' are the bosonic phonon operators and (? Pz , q is the 
phonon- matrix element. The notation A indicates that 
the factorized single-particle contributions (subscript 
S) are removed, e.g. 
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The two-particle correlations not only influence the pop- 
ulation dynamics, they also determine T k which leads to 
the decay of P. 

If one ignores r k , the optically generated state 
l^coh(O) can easily be obtained from the well-known 
coherent-limit results Q of Eqs. CD - @- For ex- 
citations resonant with the excitonic state <p v {k), we 
find P k oc <t> v (k) ^ k , 0„(k') and the conservation law 
(n k - \f + |P k | 2 = \ where n k = = /£. These co- 
herent excitons have a strict s-like symmetry since only 
then J2 k , (k') is nonvanishing. In addition, |\Pcoh(*)) 
is a Slater determinant \^ C oh(t)) = Ilk -^1(01*0) = 
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D(X) Yik at kl^o) where |^o) i s the unexcited semicon- 
ductor while Ll(t) = e i¥,k (*W/3 k (» o| k + cos/3 k (i) a^ k 
is a fermion operator with sin 2 /3 k (i) = n k (i) and e^* 1 -*' = 

Mt)/\Mt)\- 

For our subsequent discussions, it is conve- 
nient to introduce an exciton operator X„, q = 
^^(kja^^^a^k+q, with a center-of-mass mo- 
mentum q and q e (^) = m e (? l )/(m e + m/j)q. By choosing 
X = X u . and 0„(k) = /3 k e^ k , |* coh (i)) follows from 
D(X) = e x ^~* acting on the full valence band. Since 
D(X) is formally analogous to the coherent state gener- 
ator of bosonic fields we may interpret \^> C oh(t)) as 
a coherent exciton state even though X is not bosonic 

To study, how efficiently |^ C oh(t)) can be converted 
into incoherent excitons, we need to solve the full 
Eqs. CJ-©. For resonant excitation, both c Vt x,\,c and 
P k,q convert a coherent excitonic polarization into inco- 
herent two-particle populations where c qk k = c q ' k c ' k 
describes incoherent excitons. Using the cluster expan- 
sion pJ|. we derive equations for the two-particle corre- 
lations: 
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where e qk ,k is the renormalized kinetic energy of the 
two-particle state and S contains Coulomb induced in- 
and-out scattering of single-particle quantities. The 
Coulomb sums with the phase-space filling factor (f — 
f e — f h ) describe the attractive interaction between elec- 
trons and holes, allowing them to become truly bound 
electron-hole pairs, i.e. incoherent excitons which can be 
probed via the THz induced j ■ A coupling. The G term 
contains the same c V: \^\ :C and T> v ' c correlations as T in 
Eq. J3J), showing how coherent excitons are converted 
into incoherent ones. The remaining two-particle con- 
tributions are denoted as -D ros t while T symbolizes the 
three-particle Coulomb and phonon terms treated here 
at the scattering level. This way, we fully include one- 
and two-particle correlations and obtain a closed set of 
equations providing a consistent description of optical as 
well as THz excitations in semiconductors. 



The polarization to population conversion effi- 
ciency is determined from the density of incoher- 
ent ^-excitons [l2J An v — ^ q A{Xl^X u ^ = 

S k k ' *K/ (k)<^i/(k)ey _q '" +qe where C d is the normal- 
ization volume. We evaluate the conversion efficiency by 
numerically integrating the complete set of equation for a 
planar arrangement of identical quantum wires. Later on, 
we estimate the conversion efficiency also for a quantum 
well. We choose standard GaAs-type parameters and the 
wire and well sizes are taken such that the energy separa- 
tion between the two lowest exciton states is 5 meV. The 
lattice temperature is assumed to be I K such that it is 
sufficient to include only acoustic phonons 13J . To study 
the generation of incoherent excitons in their different 
quantum states, we assume 4 ps pulsed optical excitation 
resonant with either the Is- or 2s-resonance. We repeat 
the computations for different pump intensities and eval- 
uate the final quasi-stationary exciton fraction An^/n 
relative to the generated carrier density n = -i^ ^2 k / k . 
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FIG. 1: (a) For excitation at the Is exciton resonance with 
a 4 ps laser pulse (dot-dashed line), the temporal evolution 
of the induced optical polarization \P\ 2 (shaded area), to- 
gether with the generated incoherent Is (dashed line) and 2p 
(solid line) exciton densities [10 4 cm -1 ] are shown. The inset 
shows the pump (shaded area) and linear absorption (solid 
line) spectra; Ei a is the ls-exciton energy, (b) The polariza- 
tion to population conversion efficiency for Is (dashed line) 
and 2p excitons (solid line) is plotted as function of excitation 
density n. The arrow indicates the density at which the dy- 
namics is shown in a). The shaded area represents the result 
obtained without the phonon scattering. 

In Figs.nand|21 we present numerical results for pump- 
ing at the Is- and 2s-resonances, respectively. The insets 
show the spectral excitation conditions. In Fig. ^a), we 
plot the temporal evolution of the pump pulse, of the in- 
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duced optical polarization, and of the generated Is and 
2p exciton density. Figure ^b) presents the relative per- 
centage of excitons in the different quantum states show- 
ing, not surprisingly, that for Is excitation the optical 
polarization is mainly converted into incoherent Is exci- 
tons; Ani s /n is well above 90 % for low densities. This 
large conversion fraction is expected since coherent and 
incoherent Is excitons have an excellent energetic match. 
However, the generated exciton population drops well be- 
low 40 % already at moderate densities above 10 5 cm -1 
where the phase space filling factor (1 — f e — f h ) peaks 
around 0.5. For even higher densities, An^/jic/, vanishes 
since excitons start to ionize. A computation, where 
phonon scattering V is omitted, indicates that Coulomb 
scattering alone would lead only to roughly 15 % exciton 
population. 
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FIG. 2: Same as Fig .1 but for excitation at the 2s exciton 
resonance (inset), (a) Dynamics of optical polarization \P\ 2 
(dot-dashed line) and incoherent densities of 2s (shaded area) 
and 2p (solid line) excitons [10 4 cm' 1 ], (b) Conversion effi- 
ciency for Is (dashed line), 2p (solid line), and 2s (shaded 
area) excitons as function of excitation density n. 



The 2s-excitation results are presented in Fig. where 
Fig- 13b) shows An2 S /n, An2 P /n, and Ani s /n. For not 
too high excitation densities, we observe that the 2s- 
polarization is converted into a mix of 2s and 2p popu- 
lations. Whereas the amount of 2s population decreases 
monotonously with increasing excitation, the 2p popula- 
tion first increases up to 40 % before it also decreases at 
higher densities where formation of ls-excitons gradually 
becomes relevant. 

Before we analyze the physical mechanism responsi- 
ble for the significant formation of a 2p-exciton popu- 
lation, we study its signatures assuming a THz probe. 
For this purpose we evaluate the j ■ A terms in Eqs. 



and Q to compute the generated THz current Jthz = 
\3>S^)fk w ith both coherent and incoherent con- 
tributions [lj]. We determine the linear THz gain g(u>) — 
— Im [JtHz(w)/ assuming a 150 fs THz probe 

pulse capable of resolving temporal snapshots of g(co) 
during the exciton formation process. Figure |3Ji, shows 
that for ls-pumping the corresponding g{uj) is always 
negative, i.e. we find THz absorption peaked around the 
ls-2p transition. These results, as well as the asymmetric 
line shape due to transitions from the Is to energetically 
higher bound and unbound states agree well with exper- 
imental findings reported in Ref. [5j. 
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FIG. 3: Terahertz gain spectra g(uj) for different time de- 
lays corresponding to the conditions of (a) Is excitation as 
in Fig. [Ql and (b) 2s excitation as in Fig. [5^. All curves are 
identically scaled but shifted with respect to another. Here, 
E21 is the energy difference between the Is and 2p exciton 
states. 

The same analysis is repeated for the 2s excitation 
(Fig. [3d) showing that g(uj) rapidly changes from absorp- 
tion to gain. At very early times, the system consists 
mainly of coherent 2s excitons such that only absorp- 
tive transitions from the 2s to higher states are present. 
However, as incoherent 2p excitons are generated, pro- 
nounced THz gain develops at the 2p to Is transition as 
a consequence of the population inversion between these 
states. 

Now, after we have presented examples of the fully 
numerical evaluation of our many-body theory, we want 
to analyze the relevant physical mechanisms. At first 
sight the significant generation of p-type excitons might 
be unexpected since it involves a symmetry change of the 
optically generated s-type polarization. To identify the 
microscopic origin of this process, we numerically study 
the relative importance of different contributions to the 
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full theory. We find that for 2s pumping a switch off of 
phonon-induced scattering leaves the generated exciton 
fraction practically unchanged. Hence, in contrast to Is 
pumping, the conversion of a 2s polarization predomi- 
nantly results from Coulomb scattering. 

In order to understand how the Coulomb interaction 
induces symmetry changes in the polarization conversion 
we now investigate the scattering (r) and conversion (G) 
mechanisms which stem from the same fermionic corre- 
lation c v ^\,x.c — ^(o-t, a x a ^ a c) between polarization and 
fermionic density. 

(1) In a first step, we make the lowest level approxima- 
tion by replacing the microscopic Tk by a phenomeno- 
logical decay: I\ = — VP^. In this case, G q ' k ,k = 
+2P*rP k < 6q t0 and [P£Pv + E q c x k ' k ] is a constant 
of motion with respect to the scattering. Therefore, this 
simple dephasing model only converts s-type polarization 
to s-like exciton populations, in contrast to our micro- 
scopic results. 

(2) Looking at the process of excitation induced de- 
pahsing we see that Coulomb induced dephasing is 
actually a diffusive redistribution of the microscopic po- 
larizations since Ek I\ = and Ek k , G qk k = 0. At 
the same time, Eqs. © and (0 impose a strict micro- 
scopic connection between T and G. In order to analyze 
the consequences of these fundamental restrictions, we 
use a somewhat reduced model by utilizing a simplified 
structural form of the second-Born solution of Cu,a,A,c 

via Ea c^[t = -(Ex < M tt = iF(f, i>)k<, q (K- q - 
P q ) where F(f, P) is a real- valued, nonlinear functional 
of / and P containing the energy and momentum con- 
servation aspects of the Coulomb scattering. The re- 
duced model implies T k cd = — E q ^q(Pk — flc-q) an d 

G? e 'd'' k = (*K - JE-jCqfiV+q - fl*) , where C/ q = 
2Vq Ek' F(f, P)k',q- As a result, r k d removes polar- 
ization from the state Pk and redistributes it to Pk- q 
while the rate of conversion to the exciton state v be- 
comes G^^ v — |M V) q| 2 t7q. Here, the scattering ma- 
trix element is M v ^ = Ek C( k ) [Pk+ q „ - Pk- q J in- 
dicating that Coulomb scattering leads to the genera- 
tion of excitons with finite momenta whereas no pop- 
ulation in the q = state is produced. For low to 
moderate 2s excitation, we may use the approxima- 
tion Pk oc 02s (k). With the help of the symmetries 
2s (-k) = 02s(k) and 2p (-k) = -0 2p (k), we find 
M 2piq oc Ek a </ , 2 P (k)</'2s(k-|-qA) which is clearly nonzero 
for q ^ 0. 

For 2s pumping, the energy conservation aspects of Uq 
are practically the same for 2s and 2p since these states 
are nearly degenerate. As a result, the overlap of the 
wavefunctions with shifted arguments in M„ fl determines 
the conversion rate such that |Mi, iq | 2 can be used to esti- 
mate the ratio of generated 2s and 2p populations. Using 
the low-density exciton wavefunctions, we obtain a ratio 
of 1.36 of 2s over 2p population for the quantum wire, 



which is close to the numerical result in Fig. 03. Repeat- 
ing the same calculations for the quantum well, we get 
a ratio 0.99 showing that the generation of p-like states 
is strong and qualitatively similar for quantum wells and 
wires. 

Since the Coulomb interaction conserves the angu- 
lar momentum, one may ask how this conservation law 
is fulfilled when a 2s polarization is converted into 2p 
excitons. Without THz fields, Eq. (@J implies corre- 
lations with a functional dependence c q,k ~^M+ie = 
c(|q|, |k'|, |k|,cos<£k,q,cos¥>k', q ) where ¥>k(k'),q is the an- 
gle between k (k') and q. Consequently, the gener- 
ated excitons have an angular dependency cos m (^k.q = 
(e+^k.c, + e - lVk ,c,)m/2 m such that the eigenfunctions 
e +tmip an( j e -vnup Q f ^ e Tj z state are generated with equal 
probability. As a result, the total (Li) always vanishes 
such that the total angular momentum is fully conserved 
even when 2p excitons are generated. 

In summary, our microscopic study predicts significant 
formation of excitons in 2p states for excitation at the 2s 
resonance of the absorption spectrum. As a consequence, 
an exciton population inversion between the 2p and Is 
states is realized leading to gain for the corresponding 
THz frequency. Using different semiconductor materials 
this scheme may become useful for THz amplification in 
a wide spectral range. 
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